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Collisional properties of cold spin-polarized metastable neon atoms 
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We measure the rates of elastic and inelastic two-body collisions of cold spin-polarized neon 
atoms in the metastable 3 P2 state for 20 Ne and 22 Ne in a magnetic trap. From particle loss, we 
determine the loss parameter of inelastic collisions (3 = 6.5(18) x 10 -12 cm 3 s -1 for 20 Ne and /3 = 
1.2(3) x 10~ n cm 3 s _1 for 22 Ne. These losses are caused by ionizing (i.e. Penning) collisions and occur 
less frequently than for unpolarized atoms. This proves the suppression of Penning ionization due 
to spin-polarization. From cross-dimensional relaxation measurements, we obtain elastic scattering 
lengths of a = -180(40) a for 20 Ne and a = +150±|g a for 22 Ne, where a = 0.0529 nm. 

PACS numbers: 34.50.-s, 34.50.Fa, 32.80.Pj, 82.20.Pm 



The experimental investigation of the collisional inter- 
action of cold and ultracold atoms has been a driving 
force for important new physics in recent years. The 
successful demonstration of Bose-Einstein condensation 
(BEC) in dilute atomic gases certainly is one of the most 
prominent examples of this development 0. Other ex- 
amples are the realization of Fermi-degenerate systems 
(see for an overview) or the generation of molec- 
ular BECs 3J. The collisional properties of ultracold 
atoms critically determine the achievability of quantum- 
degenerate systems, their stability, their coherence prop- 
erties, and the characteristics of collective phenomena 
(see jj, Li( for a review) and thus deserve detailed inves- 
tigation. All along, experimental work on cold collisions 
for an additional atomic species has initiated a new line of 
exciting research. Special interest nowadays arises from 
the investigation of non- alkali systems, such as He* 
Yb 0, Cr 8, 9], with their specific interaction properties. 

An important novel class of atomic species in this re- 
spect are cold and ultracold rare gas atoms in metastable 
triplet states (RG* atoms). Due to their high internal 
energy, these atoms show collisional properties and ex- 
perimental provisions absent in other atomic systems. 
Important results have already been achieved in He*, 
such as the realization of BEC |J 0, the implementa- 
tion of highly efficient electronic detection of quantum- 
degenerate atomic samples its application to the de- 
termination of the elastic scattering length ^(| > an d the 
formation of He*, dimers in purely long-range states . 
With this paper, we extend the experimental investiga- 
tion of cold collisions of metastable atoms to the case of 



neon. 



Compared to helium, neon has a more complex inter- 
nal structure. This causes for example the occurrence of 
more than one metastable triplet state and of anisotropic 
electrostatic interaction potentials, which influence the 
dynamics of elastic and inelastic collisions 0, 0] . In- 
elastic interactions of RG* atoms are particularly inter- 
esting due to their specific exoergic character. With high 
probability, RG* atoms exhibit Penning ionization reac- 



tions 14], where one of the atoms is deexcited to the 
ground state, the other atom is ionized, and an elec- 
tron is released. Penning ionization may be suppressed if 
the colliding atoms are spin-polarized to a spin-stretched 
state since then Penning ionization does not conserve the 
spin quantum number. The extent of suppression de- 
pends critically on anisotropic contributions to the inter- 
action. In the case of He* these contributions are small 
and ionization is suppressed by four to five orders of mag- 
nitude [ 101 ] . The heavier rare gases however have stronger 
anisotropic contributions, since the excitation of an elec- 
tron to the metastable state creates a p 5 -core with non- 
zero orbital angular momentum. Indeed, no suppression 
of Penning ionization was observed in Xe* 0] . For neon, 
Doery et al. calculate suppression ratios between 1 and 
10 4 , which depend most sensitively on exact interaction 
potentials 01 > which were not available at the time. Ex- 
perimental investigations are needed, but so far measure- 
ments of the ionization rate have only been reported for 
unpolarized atoms [lfl. UtI Il8j | . 

In this paper, we present the experimental determina- 
tion of the rates of elastic and inelastic collision of spin- 
polarized bosonic 20 Ne and 22 Ne atoms in the J=mj=+2 
substate of the metastable 3 P2 state. For both isotopes, 
we determine the elastic scattering length a, the loss pa- 
rameter of inelastic collisions (3, and demonstrate the sup- 
pression of Penning ionization due to spin-polarization. 

In our experiment, we capture Ne* atoms in a 
magneto-optical trap from a Zeeman-decelerated atomic 
beam 0. After optical pumping into the mj= + 2 state, 
the atoms are transferred into a magnetic Ioffe-Pritchard 
trap with radial gradient of 295G/cm and axial curva- 
ture of 215G/cm . As a last step of preparation, we 
apply one-dimensional Doppler cooling [20j, where the 
trapped atoms are irradiated by two a + polarized laser 
beams along the symmetry axis of the magnetic trap. 
For this purpose, we operate the trap with an offset mag- 
netic field of 25 G. The vibrational frequencies, measured 
with 20 Ne, are u> x = 2ir X 80(1) s _1 in axial direction and 
u) r = 2tt x 186(1) s -1 in radial direction. By this prepa- 
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FIG. 1: Cross-dimensional relaxation of the aspect ratio of 
22 Ne ensembles for different initial densities n(0). Higher 
densities lead to shorter relaxation times 7"*. The lines are 
fits according to Eq. During relaxation n and T change. 
We compensate for the resulting change in 7 rc i by using the 
rescaled time t*(t) = f* (n(t') v(t')) / (n(0) v(0)) dt', similar 
to [ae- 
ration we obtain typical temperatures of T x = 450 /iK 
in axial and T r — 600 /iK in radial direction. A typical 
atom cloud contains N = 2 x 10 s atoms, at a mean den- 
sity n — 4 x 10 10 cm -3 , the latter being the ensemble 
average of the number density n(r): n — N^ 1 Jn 2 (r) d 3 r. 
By a Stern-Gerlach experiment, we could verify that the 
spin-polarization persists during Doppler cooling. We 
could not detect any atoms in the J=2, mj<l substates, 
from which we deduce a lower bound of 95% for the atom 
population in the m,/=2 substate. 

As one-dimensional Doppler cooling in the magnetic 
trap gives different temperatures, T x < T r , we obtain 
an ensemble out of equilibrium. In each dimension how- 
ever, potential and kinetic energy are balanced such that 
oj x r r \ cr x (r) = ^B^^/m, where u x ( r ) is the spatial 
width in axial (radial) direction and ^/feeT^^/m the 
velocity spread, with Boltzmann's constant fee, and the 
atomic mass m. After cooling, the trapped ensemble 
trends back to equilibrium, which can be observed by 
a changing aspect ratio of the spatial distribution. In 
parallel, the number of atoms decreases with time and 
their mean temperature increases. We observe heating 
rates which are density dependent and range from a few 
/iK/s up to 50/iK/s for increasing densities. 

In general, collision rates depend on the atom cloud's 
density and temperature. We derive n, T x , and T r from 
absorption images which are taken between 0.5 and 2 ms 
after switching off the magnetic trap. For the determi- 
nation of these parameters, statistical uncertainties (typ- 
ically 2 — 3%) are negligible as compared to systematic 
uncertainties. Uncertainties are given as one standard de- 
viation and represent the quadrature sums of statistical 
and systematic contributions, dominated by the system- 
atic uncertainty in density (23% < An/n < 28%). 

We determine the rates of elastic collisions from cross- 
dimensional relaxation measurements 0, 0, |2^]. The 
equilibrating transfer of energy from the radial to the 



FIG. 2: Effective relaxation cross sections <j rc i obtained 
from cross-dimensional relaxation measurements and calcu- 
lations: (left) For 20 Ne we obtain a scattering length of 
a — —180(40) ao. The measurements between 500 and 700 /iK 
are also consistent with a positive a — +22(5) ao; (right) The 
large <7 rc i-values of 22 Ne are close to the unitarity limit. Our 
data agree best with theory for a — +I50I5Q a o- We can rule 
out negative scattering lengths since theory cannot match our 
data, as exemplified for a = — 6500 ao (see text). 

axial dimension changes the atom cloud's aspect ratio 
A(t), according to 

A(t) = - 7 rcl(A(<) - A cq ) . (1) 

Here, A eq is the equilibrium aspect ratio and 7 re i the re- 
laxation rate. We observe reequilibration of 20 Ne and 
22 Ne (Fig. nj for different initial densities n(0) and find 
relaxation rates which are directly proportional to n(0) 
[l8| . Hence relaxation is predominantly caused by col- 
lisions and ergodic mixing is of no relevance for our de- 
termination of 7 ro i 23]. Therefore, we can describe the 
relaxation rate of the atomic ensemble by an effective 
relaxation cross section cr re i , 

7rd = a- rc \nv , (2) 

where v = (l6fcsT /urn) is the average relative veloc- 
ity of the colliding atoms, with the mean temperature 
T = i{T x + 2T r ). 

Since the relaxation process is driven stronger by col- 
lisions with high relative velocity, Kavoulakis et al. have 
derived a rc \ oc (a e i(v ) u) n st. from a non-standard ther- 
mal average of the elastic collision cross section a e \(v), 
(Eqs. (72) and (88) in 24]), where v denotes the rel- 
ative velocity of the colliding atoms. Thereby cr re ] be- 
comes a function of scattering length a and temperature: 
Crci = CT ro i(a, T). 

Figure [21 shows the relaxation cross sections of 20 Ne 
and 22 Ne as a function of temperature [25i |. We have 
prepared atom clouds at different temperatures either by 
varying the efficiency of Doppler cooling, or by an ad- 
ditional adiabatic expansion within the trap. The tem- 
perature assigned to each measurement is the average of 
temperatures T(t) weighted with A(t). This procedure 
reflects that the period of fastest ch ang e in A(t) is most 
important for the extraction of 7 re i [26j . 
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In our measurements, relaxation is dominated by s- 
wave collisions between spin-polarized atoms, as the cen- 
trifugal barrier for d- waves is fee x 5.6 mK in energy. In 
the temperature range we explore, the elastic scattering 
cross section depends on both the magnitude and the 
sign of the scattering length a, as opposed to the case of 
ultracold collisions, where a c \ — > 87ra 2 . To analyze the 
data, we use the temperature dependence of <j ro \(a,T) 
which we derive from numerical calculations: we start by 
calculating a c i(v) as a function of v for different input 
values of a by solving the radial Schrodinger equation. 
With the non-standard average {cr c i(v) v) ns t. we calcu- 
late cr re i(a, T) for temperatures between 100 and 800 /xK 
(lines in Fig. EJ and determine the scattering length by 
fitting these curves to the experimental data (x 2 mini- 
mization procedure). 

These calculations are based on the short-r ang e ab- 
initio potentials given by Kotochigova et al. [221], and 
the long-range van-der-Waals potentials given by Dere- 
vianko et al. We treat the collision process by a 

single-channel calculation, instead of including the five 
potentials involved, for three reasons: (i) the sensitivity 
of the u-dependence of a e Av)on the short-range poten- 
tials is negligible (see also |8j, |lj| ) ; (ii) the 4% uncertainty 
in the van-der-Waals coefficients for the different poten- 
tials involved is comparable to their differences |28| : (iii) 
the variation in the calculated <7 re i(a, T) caused by these 
differences are negligible when compared to experimental 
uncertainties in the measurement of a rc \. 

For 20 Ne, we obtain a negative scattering length of 
a = -180(40) ao, with Bohr's radius a = 0.0529nm. 
Between 500 and 650 /iK, our data are also consistent 
with a positive value of a = +22(5) ao- The sign of the 
scattering length is determined by the measurements at 
T w 200 fiK, giving a 74 times higher probability for 
a = — 180 ao than for a — +22 ao- For 22 Ne, we obtain 
relaxation cross sections cr ro i which are three to five times 
larger as compared to 20 Ne (Fig. EJ. We obtain a — 
+150^50 ao. Due to the proximity of the cr rc i-data to the 
unitarity limit, the extracted a is less precise than in the 
case of 20 Ne. Our calculations show that the relaxation 
cross sections for a < do not match our data. We 
exemplify this by including a plot of <r re i for a = —6500 ao 
in Fig. |21 All curves calculated for —6500 a < a < lie 
below this curve, and thus disagree with our data. We 
therefore conclude that the scattering length of 22 Ne is 
positive. 

From atom loss in the magnetic trap, we determine 
the rate of binary inelastic collisions of a spin-polarized 
ensemble. The decrease in atom number N(t) (Fig. El 
left) stems from one-body losses (oc N) and two-body 
losses (oc TV 2 ). We have no signature of higher order loss 
processes in our measurements. One body-losses with 
rate a are mainly caused by the 14.73 s lifetime of the 3 P2 
state 19] and by background gas collisions. Two-body 
losses depend on the probability of two atoms to collide 
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FIG. 3: (left) Atom loss of a spin-polarized 20 Ne ensemble. 
Initially, 2 x 10 s atoms were trapped at a mean density n = 
2 x 10 10 cm~ 3 . (right) Data presented according to Eq. Q 
and a linear fit giving a one body loss rate a = (10.3(3) s) -1 
and a loss constant — 6.5(18) x 10" 



(oc n 2 (r)) and on the loss parameter (3. The decay of the 
particle number is described by 



dt 



N(t) = -aN(t) - (3 



N 2 {t) 
V cS (t) 



(3) 



with the effective volume V e g = N/n. We find that V e s 
increases due to heating, so that an explicit solution of 
Eq. J2J requires assumptions on the heating mechanisms 
involved. We circumvent this by formally integrating 
Eq. Js3 an d obtain a linear relation between quantities, 
which are easily derived from measured data: 



1 N(t) 



n(t')dt' 



(4) 



The constant offset in this relation is the one-body loss 
rate a and the slope is the two-body loss parameter (3. 

By presenting the data of a trap loss measurement 
(Fig. El (left)) according to Eq. |@J|, we find the expected 
linear dependence (Fig. El (right)). From this, we get the 
one-body decay rate a -1 = 10.3(3) s, and the two-body 
loss parameter (3 = 6.5(18) x 10~ 12 cm 3 s _1 . The uncer- 
tainty in (3 is dominated by the systematic uncertainty 
in n. Comparing (3 to the loss parameter of an unpo- 
larized ensemble /3 u „ P oi = 2.5(8) x lO-^cnrV 1 [Hill, 
we obtain the suppression of inelastic collisions due to 
spin-polarization f3 unpo \/ f3 = 38(16) for 20 Ne. 

We get a second independent measurement of the rate 
of inelastic collisions from observed heating rates. Be- 
sides other possible heating mechanisms, two-body loss 
leads to intrinsic heating: these collisions happen most 
likely in the center of the trap. Therefore, the average 
energy carried away by a lost atom is less than the mean 
energy per trapped atom. By fitting the resulting heat- 
ing rate T/T — (3n/4 29] to the rising mean temper- 
ature T, we obtain /3 hc at = 7(2) x 10~ 12 cm 3 s _1 . This 
value agrees well with the result from the trap loss mea- 
surement and indicates that other heating mechanisms 
do not contribute significantly to observed heating rates. 
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TABLE I: Summary of collision parameters: effective relax- 
ation cross section a ro i for T = 200 and 550 /iK, scattering 
length a, two-body loss parameters f3 and /3 unp oi for polarized 
and unpolarized atoms and their suppression ratio. 



J Ne 



2 Ne 



<Trcl(~ 200^iK) 

(« 550^K) 

a 

P 

/^unpol 
/3unpol/ f3 



[io- 
[io- 



[10~ 12 cmV 
[10 -12 cmV 



8(2) 
2.8(7) 
-180(40) 
6.5(18) 
250(80) 
38(16) 



30(8) 
13(3) 
+150±«g 
12(3) 
80(50) 
7(5) 



We performed similar measurements with 22 Ne at an 
initial atom number of 9 xlO 7 (n = 1.4 x 10 10 cm -3 ) yield- 
ing a one-body decay rate of a -1 = 14.5(2) s. We find 
a loss parameter of /3 = 1.2(3) x 10~ n cm 3 s~ 1 from the 
decay measurement, and /3 nC at = 1-3(3) x 10~ n cm 3 s _1 
from observed heating. With a loss parameter (3 unpo i — 
8(5) x 10 -11 cm 3 s _1 [lg for unpolarized atoms, the sup- 
pression of two-body losses is /3 unpo \/f3 = 7(5) for 22 Ne. 

In an additional experiment, we simultaneously de- 
tected atoms with absorption imaging and measured the 
rate of ions escaping from a cloud of trapped 22 Ne atoms 
with a multichannel plate (MCP) detector. From these 
measurements (analogous to |l9j ) we can establish the 
ratio of the rate of ionizing inelastic collisions (/3i on n) to 
the total rate of inelastic collisions (/3n): Aon//? = 1-1(2), 
where the uncertainty is dominated by the calibration 
uncertainty of the MCP detector. These measurements 
confirm, that the observed inelastic collisions are caused 
by Penning ionization to a very high degree (> 90%). 

Our measurements prove the suppression of Penning 
ionization due to spin-polarization. Compared to He* 
this suppression is much less pronounced. The observed 
values for /3 U npoi//3 do not confirm that large suppres- 
sion ratios are likely [12T l30| . We expect that refined 
calculations of the ionization rate with precise interac- 
tion potentials |27| 28J will reproduce the experimental 
results. 

As a summary, Table H] gives the collisional parame- 
ters for both neon isotopes as presented in this paper. 
With these measurements, a significant contribution to- 
wards the understanding of the collisional properties of 
spin-polarized metastable neon atoms is gained. These 
results will initiate further theoretical and experimental 
investigations, such as detailed calculations of the molec- 
ular potentials and their manipulation by external elec- 
tromagnetic fields, e.g. to modify the rate of elastic or 
inelastic collisions. Concerning the quest for BEC, the 
determined ratio of elastic to inelastic collisions suggests 
22 Ne to be better suited for evaporative cooling than 
20 Ne. Accordingly, we have carried out first experiments 
of evaporative cooling of 22 Ne and an increase in phase 
space density was already observed 01 • Whether a BEC 



of 22 Ne can be realized hereby, will be revealed by forth- 
coming evaporation experiments. 
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